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Abstract 

We observe the dynamic formation of Cs2 molecules near Feshbach resonances in a cold sample 
of atomic cesium using an external probe beam. This method is 300 times more sensitive than 
previous atomic collision rate methods, and allows us to detect more than 20 weakly-coupled 
molecular states, with collisional formation cross sections as small as a = 3 x 10~ 16 cm 2 . We 
propose a model to describe the atom- molecule coupling, and estimate that more than 2 x 10 5 Cs2 
molecules coexist in dynamical equilibrium with 10 s Cs atoms in our trap for several seconds. 

PACS numbers: 34.50.-s, 05.30. Jp, 32.80.Pj, 67.40.Hf 
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The newly emerging field of cold molecules provides intriguing possibilities in ultracold 
gas studies. In the realm of "super-chemistry" , phase-coherent chemical reactions between 
particles in well-defined internal and external states are dominated by quantum statistics 
and many-body effects PJ. The coherent nature of the atom-molecule coupling has been 
recently observed in Bose-Einstein condensates 0. Quantum manipulation of molecular 
states has also been proposed for quantum computation ||. Finally, in precision tests of 
time-reversal symmetry based on matter-wave interferometry, polar molecules promise a 
sensitivity several orders of magnitude higher than atoms |4[j. 

However, the production of cold molecules presents a major experimental challenge. Di- 
rect laser cooling of molecules has been proposed using either multiple laser frequencies @ 
or an external cavity |J; however, these methods involve exquisite experimental control, and 
have not yet been implemented. Photoassociation of cold atoms into molecules is limited by 
the phase space density of the initial atomic sample 0, and consequently only a few thou- 
sand molecules have been produced from thermal samples with this method. Two-photon 
photoassociation from an atomic Bose-Einstein condensate has been limited to similarly 
small numbers of molecules due to large atom-molecule inelastic collision losses 0. Alter- 
native routes to cold molecules by slowing down a supersonic molecular jet ||, or by buffer 
gas loading of a magnetic trap [JHJ , can collect large numbers of molecules, but the resulting 
temperatures are much higher than in the approaches using cold atoms. 

In this work, we employ a radiation-induced loss mechanism to monitor the dynamic 
formation of cold Cs2 molecules near Feshbach resonances. Our method is similar to the 
detection technique introduced by Heinzen and coworkers however, instead of photoas- 
sociating ground state atom pairs, which requires knowledge of the excited-state molecular 
level structure, we use a blue-detuned probe beam to photodissociate them. In a dense sam- 
ple of optically trapped cesium atoms we observe a class of weakly-coupled Cs 2 molecular 
ground states which, in the absence of the probe light, are collisionally metastable at the 
density of our samples. Since these states are long-lived, a large molecular population can 
build up even in a thermal sample. 

Our method for detecting the formation of ground state molecules, shown in Fig. |1], can 
be understood as a radiative collision process [12, For a typical detuning of A/2tc ~ +1 
THz from the free-atom transition, radiative loss occurs when the dipole-dipole interaction 
Vd{r) ~ C3/V 3 between two colliding atoms shifts the excited state into resonance with the 
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FIG. 1: Illustration of radiative Feshbach spectroscopy. In the vicinity of a Feshbach resonance 
the scattering state ^f s couples to a ground molecular state \& m . A probe beam tuned to the blue 
of the free-atom transition excites the molecules, which then dissociate and are lost from the trap 
due to the energy imparted by the strongly repulsive excited-state potential. 
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laser field at the Condon point r c , where r c « 2nm satisfies HA = Vd(r c ) 
diative coupling amplitude is therefore proportional to the two-atom wavefunction at the 
Condon point, \I/(r c ). On a Feshbach resonance JH], the scattering state ^ s (r) couples 
to a molecular state *S? m (r), that has a much larger amplitude at short range and conse- 
quently at the Condon point. The coupling amplitude is then enhanced by a factor of 
A ~ |\I/ m (r c )|/| , I' s (r c )| ~ (r a /r m ) 3 / 2 , where r a (r m ) is the mean interatomic separation in 
the scattering (molecular) state. Typically, r a ~ 1/xm at our atomic density of 10 13 cm -3 . 
r m can be estimated from the ground state van der Waals interaction V g = —Cg/r 6 , and the 
binding energy of the molecular state, which is on the order of the atomic hyperfine splitting 

(C 6 /H h )V 6 ~ 2nm. 



H h |J. Given C 6 = 6859 a.u. |J and H h /h = 10GHz, we find r n 
The radiative coupling amplitude of the molecule by the far blue-detuned beam is therefore 
a factor of A ~ 10 4 stronger than that of two trapped atoms. This indicates that molecules 
can be detected efficiently without appreciably perturbing the atoms. 

Compared to lighter alkalis, cesium atoms display unusually rich cold-collision phenom- 
ena, due primarily to the large dipolar interaction HTE|| . The angular dependence of this 
interaction allows an incident partial wave with angular momentum L to couple to molecu- 
lar states with orbital angular momentum N = L ± 2nh, in nth order [|T^] . Consequently, 
many molecular states become available to scattering atom pairs when this interaction is 
strong. In the following, the above process is referred to as an L — > N resonance, e.g. an 
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FIG. 2: Radiative loss coefficient as a function of magnetic field for a mean density of 1.5 x 10 13 cm -3 , 
95% population in \F = 3, uif = 3), and a temperature of 5.5/iK. The wavelength and the mean 
intensity of the probe beam are 846.52nm and 33W/cm 2 , respectively. In the inset, the radiative 
spectrum (lower curve) shows higher sensitivity than the collision rate measurement (upper curve). 
The stars indicate |3,3) + |3,2) collision resonances. 

S — > G resonance occurs when an s-wave scattering state couples to a molecular state with 
an orbital angular momentum N = 4h. 

In this work, we apply radiative Feshbach spectroscopy to two types of cold collisions: 
When both scattering atoms are in the lowest hyperfine ground state \F = 3,mp = 3), 
denoted in the following by |3,3) + |3,3), the scattering state \T/ S can couple to molecular 
states ^> m with even orbital angular momentum as the relative energy of \T/ S and \T/ m is tuned 
by means of a magnetic field. We find multiple S — > D, S — > G Feshbach resonances, and one 
S — ► D shape resonance for samples polarized in |3,3) fL7j. For |3,3) + |3, 2) collisions, on 
the other hand, both even and odd partial waves are allowed, and we observe an abundance 



oiS^>D,S^G,P^P and P — > F Feshbach resonances |L7|. 

Our experimental setup has been described previously in Ref . [ 18 1 . About 3 x 10 8 atoms 
are loaded into a far-detuned dipole trap formed by a YAG laser at 1064nm using two 



pulses of Raman-sideband cooling |18], [0|. We then further optically pump so that up to 
98% of the remaining atoms are in |3,3), and by adjusting the optical pumping frequency 
prepare mixed samples as desired with up to 10% in |3,2). For all of the experiments 
described here, less than 1% (0.25%) of the atoms are found to occupy |3, 1) (all other 



4 




10 20 30 40 

magnetic field [G] 



FIG. 3: A detailed radiative loss spectrum with a resolution of 5mG for a probe beam at 844nm. 
90% (10%) of the atoms are in the |3, 3) (|3, 2)) state. The temperature and density of the sample 
are 3.5fiK and 1.8 x 10 13 cm -3 , respectively. "S" indicates a shape resonance; the stars indicate 
1 3, 3) + |3, 2) collision resonances. 



states), as has been determined by microwave spectroscopy ||20|| . Within 15ms the atoms 
thermalize to a temperature between 3fiK and 6/zK, and density near 10 13 cm~ 3 by means 
of elastic collisions; the atom number is measured by fluorescence detection and the density 
is then calculated from the measured temperatures and trap vibration frequencies. We 
adiabatically ramp the magnetic field from 0.1G to 2G in 1ms, and then to any desired 
field value up to 230G in another 1ms. No depolarization of the sample is observed during 
this procedure. A Titanium-Sapphire laser provides the far-detuned probe beam with an 
e~ 2 beam waist of 2.0mm (0.6mm) in the vertical (horizonal) direction, which uniformly 



illuminates the sample whose vertical (horizontal) size is 580/xm (60/im) [jl3 1 . The typical 
mean intensity and wavelength, optimized for radiative detection sensitivity, are 20W/cm 2 
and 847nm, respectively. Loss coefficients are extracted from the evolution of the atomic 
density, assuming a Gaussian distribution, according to n = —n/T — Gn 2 — (4/3) 3 / 2 i\n 3 , 
where n is the mean atomic density, T = 3 s is the one-body lifetime, and G(K) is the two 
(three)-body decay coefficient. 

Fig. |2| shows the radiative loss rate from 0G to 230G, measured with a resolution of 
lOOmG. All seven observed strong resonances are due to |3,3) + |3,3) collisions. The lower 
curve in the inset shows a higher-resolution scan performed on a sample with a larger 
impurity in |3,2). Seven additional smaller resonances are evident, which originate from 
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1 3, 3) + 1 3, 2) collisions. The upper curve in the inset is the corresponding measurement of 
elastic collision-induced evaporative loss performed under similar conditions [ p0| . 

It is worth noting that only the strongest resonance at 48.0 G could be detected with 
the collision rate measurements in Ref. [EDI. There are two reasons for this: first, there 
is no appreciable trap loss due to inelastic collisions near these resonances. In particular, 
for 1 3, 3) + 1 3 , 3) collisions, where binary inelastic collisions are endothermic and therefore 
forbidden, the resonant three-body loss is observed to be below our detection resolution of 
K 3 = 3 x 10 -27 cm 6 /s even on the strongest S — > G wave resonance at 53 G. Second, these 
resonances are sufficiently narrow in energy that the thermally averaged change in elastic 
collision rate is smaller than our detection limit of 3%. Clearly the radiative Feshbach 
resonance detection offers higher sensitivity and accuracy for locating molecular states. 

The inset to Fig. |2] shows two complementary aspects of Feshbach resonances: When the 
field shifts a molecular state into resonance with the scattering continuum, an increase in 
radiative loss indicates the appearance of a molecular population; at the same time a change 
in evaporation loss rate indicates a resonant alteration of the effective elastic collision cross 
section. The difference between these two effects is exemplified by the fact that the radiative 
resonance at 47.97 G occurs at a slightly lower field value than does the minimum in the 
evaporation loss at 48.02 G. While the radiative resonance occurs when the molecular state 
is tuned to the scattering continuum (Feshbach resonance), the minimum in the elastic 
cross-section occurs when the background scattering length is cancelled by the resonant 
contribution from the bound state (Ramsauer-Townsend effect). The finite width of the 
resonance causes these two points to be distinct. For stronger (wider) resonances this is 
more dramatic: the wide minimum in the elastic cross section observed at 17. 1G in Ref. |20[ 



is in fact induced by the Feshbach resonance at a negative field of — 8.25G (or equivalently, 
positive field for the |3, —3) state) p7| . 

We model the radiative process using a chemical rate-equation approach that treats atoms 
and molecules as two different species |H| . In a thermal sample of de Broglie wavelength X^b 
with low phase space density = n\ z dB <C 1, the evolution of the mean atomic (molecular) 
density n (m) near a Feshbach resonance can be described classically: Based on detailed 
balance between one molecular state and 1/0 scattering states in the relative coordinate, 
m = \(f>n molecules per unit volume can coexist with n atoms in thermal equilibrium in 
the absence of collision loss P2| . The formation and destruction rates of molecules are then 
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FIG. 4: Radiative loss rate at 53. 5G vs. intensity for atomic densities of n\ = 1.2 x 10 13 cm -3 (solid 
circles) and n-i = 0.6 x 10 13 cm _3 (open circles). The wavelength of the probe beam is 846. 5nm, and 
the temperature of the sample is 5.5/uK. We fit the two curves with Eq. |4|. 

described by a<pn/2 and —am, where a characterizes the strength of the Feshbach resonance. 
Given a leading order collisional loss — 7nm due to atom-molecule collisions, and radiative 
atom loss induced by the probe beam L = 25m/, where 7 and 5 are rate constants, and / is 
the probe beam intensity, we obtain the following rate equations: 



m = a <frn/2 — a m — 7 nm — L/2 (1) 
L = 25 Im. (2) 

The stationary solution for the molecular density m and the fractional atom loss rate 
L/n are then given by 



m = {(t),2)n (3) 

L /n = 5 I± (A) 

1 l + (j/a)n+(6/a)f 1 } 

To test this model, we measure the radiative atom loss as a function of probe beam 
intensity / for different atomic densities, as shown in Fig. |] for the strongest S — > G 
resonance at 53.50 G. Here the mean atomic density is switched between n\ = 1.2 x 10 13 cm~ 3 
and n 2 = 0.6 x 10 13 cm -3 by changing the trap loading parameters. The phase space densities 



are <p\ ~ 1/100 and 02 ~ 1/200, respectively. The loss constant L/n is derived from 
the measured radiation exposure time t necessary to deplete 10% of the sample, namely, 
L/n = 0.1/t. The polarization purity of 95% in |3, 3) and the temperature of 5.5/iK do not 
vary systematically by more than 5% between the two settings. We have verified that no 
additional loss is due to field ramping or radiative processes induced by the very far detuned 
dipole trap at 1064nm. 

At low probe beam intensity I < 5W/cm 2 , the fractional loss rate is linear in /. This 
suggests that the molecular and atomic populations reach a dynamic equilibrium where the 
probe beam only weakly depletes the molecules, without appreciably changing their density. 
At high probe beam intensity / > 10W/cm 2 , the fractional loss rate is proportional to 
the sample density and is independent of the probe beam intensity. This implies a quick 
depletion of molecular population by the probe beam, leading to a loss rate that is therefore 
equal to the formation rate of the molecules. As indicated by the solid lines in Fig. |4], Eq. |] 
agrees well with the experimental results, and a simultaneous fit to both sets of data with 
only two free parameters yields a = 364(20)s _1 and S/'j = 3.3(6) x 10 12 (Wcm) _1 

From the above measurements, we can first deduce the sensitivity of radiative Feshbach 
spectroscopy for detecting molecule formation. Given a, the molecule formation cross section 
at 53 G is a a = a<p/vn = 8 x 10~ 14 cm 2 , where v is the mean velocity of the atoms. While 
this resonance induces a maximum loss rate of 4s _1 , the loss rate of the weakest observed 
resonance at 7.8G is approximately 0.04s" 1 , with a signal-to-noise ratio of 3. We estimate 
that the sensitivity of our detection, expressed in terms of a molecule formation cross section, 
is then a a /300 ~ 3 x 10~ 16 cm 2 . Compared to the nominal elastic cross section for atoms in 
|3, 3) of 1.5 x 10~ 10 cm 2 , this sensitivity corresponds to the creation of less than two molecules 
every 10 6 collisions. 

Furthermore, we can calculate the stationary number of molecules in the gas from the 
measured rate constants. While the latter only reveal the formation constant a, we can 
safely put an upper limit on the collisional loss 7 from the unitarity limit of the s-wave cross 
section for inelastic collisions: a < ir/k 2 . We then estimate 7 = <j 1 v 1 < (ir/k 2 )(hk/ /i), where 
h 2 k 2 /2fi ~ \ksT is the collision energy, \i = |M the reduced mass, and M the mass of the 
Cs atom. Given Eq. |3] and a, we obtain 7 < 2.8 x 10 _11 cm 3 s _1 and a ratio of molecular to 
atomic density m/n > 2.6(5) x 10~ 3 at the atomic phase space density of <fi = 1/100 and low 
radiative excitation. Equivalently, more than 2.6(5) x 10 5 molecules coexist with 10 8 atoms 
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in our trap. 

In conclusion, we have used radiative Feshbach spectroscopy to detect more than twenty 
new Cs2 molecular states with high sensitivity and resolution. Since no strong collision 
loss accompanies these resonances, as many as 3 x 10 5 molecules accumulate and coexist 
dynamically with 10 8 atoms in our dipole trap. In the future, we hope to explore the regime 
where the atom-molecule coupling strength exceeds the thermal energy of the gas, and its 
coherent nature may become evident even in our thermal sample. 

This work was supported in parts by grants from AFOSR and the NSF. 



[1] D.J. Heinzen, Roahn Wynar, P.D. Drummond and K.V. Kheruntsyan, Phys. Rev. Lett. 84, 
5029 (2000). 

[2] E.A. Donley, N.R. Claussen, S.T. Thompson, and C.E. Wieman, |cond-mat/0204436| . 

[3] D. DeMille, Phys. Rev. Lett. 88, 067901 (2002); C. Chin, V. Vuletic, A.J. Kerman and S. Chu, 

in Proceedings of the XV. International Conference on Laser Spectroscopy (World Scientific) 

(2001). 

[4] O.P Sushkov, V.V. Flambaum, Sov. Phys. 48, 608 (1978); E.D. Commins, Adv. in At. Mo. 

and Opt. Phys. 40, 1 (1999). 
[5] H. Wadi and E. Pollak, J. Chem. Phys. 110, 11890 (1999). 
[6] V. Vuletic and S. Chu, Phys. Rev. Lett. 84, 3787 (2000). 

[7] J.D. Miller, R.A. Cline and D.J. Heinzen, Phys. Rev. Lett. 71, 2204 (1993); B.L. Tolra, C. 

Drag, and P. Pillet, Phys. Rev. A 64, 061401 (R) (2001); A. N. Nikolov, J. R. Ensher, E. E. 

Eyler, H. Wang, W. C. Stwalley, and P. L. Gould, Phys. Rev. Lett. 84, 246 (2000). 
[8] R. Wynar and R.S. Freeland and D.J. Han and C. Ryu and D.J. Heinzen, Science 11, 1016 

(2000); J.M. Gerton and D.V. Strekalov and I.D. Prodan and R.G. Hulet, OSA, Quantum 

Electronics and Laser Science Conference, 2000. 
[9] H.L. Bethlem, G. Berden and G. Meijer, Phys. Rev. Lett. 83, 1558 (1999). 
[10] J. D. Weinstein, R. deCarvalho, K. Amar, A. Boca, B. C. Odom, B. Friedrich, and J. M. 

Doyle, J. Chem. Phys. 109, 2656 (1998). 
[11] Ph. Courteille, R. S. Freeland, D. J. Heinzen, F. A. van Abeelen and B. J. Verhaar, Phys. 

Rev. Lett. 81, 69 (1998). 



9 



[12] K. Burnett and P.S. Julienne and K.-A. Suominen, Phys. Rev. Lett. 77, 1416 (1996). 

[13] V. Vuletic, C. Chin, A.J. Kerman and S. Chu, Phys. Rev. Lett. 83, 943 (1999). 

[14] E. Tiesinga and A.J. Moerdijk and B.J. Verhaar and H.T.C. Stoof, Phys. Rev. A 46, R1167 

(1992); E. Tiesinga and B.J. Verhaar and H.T.C. Stoof, Phys. Rev. A 47, 4114 (1993). 
[15] Paul J. Leo, Carl J. Williams, and Paul S. Julienne, Phys. Rev. Lett. 85, 2721 (2000). 
[16] F.H. Mies, C.J. Williams, P.S. Julienne, and M. Krauss, J. Res. Natl. Inst. Stand. Technol. 

101, 521 (1996). 
[17] Paul Leo and Eite Tiesinga, private communication. 

[18] V. Vuletic, C. Chin, A.J. Kerman and S. Chu, Phys. Rev. Lett. 81, 5768 (1998). 

[19] A.J. Kerman, V. Vuletic, C. Chin and S. Chu, Phys. Rev. Lett. 84, 439 (2000). 

[20] C. Chin, V. Vuletic, A.J. Kerman and S. Chu, Phys. Rev. Lett. 85, 2717 (2000); V. Vuletic, 

A.J. Kerman, C. Chin and S. Chu, Phys. Rev. Lett. 82, 1406 (1999). 
[21] E. Timmermans, P. Tommasini, M. Hussein and A. Kerman, Phys. Report 315, 199 (1999). 
[22] Maximization of the classical entropy yields the same result. 



10 



